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Response of fermions in Chern bands to spatially local quenches 
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We study the dynamical evolution of Chem-band systems after subjecting them to local quenches. For 
open-boundary systems, we show for half-filling that the chiral nature of edge states is manifested in the time- 
dependent chiral response to local density quenches on the edge. In the presence of power-law traps, we show 
how to mimic the half-filling situation by choosing the appropriate number of fermions depending on the trap 
size, and explore chiral responses of edges to local quenches in such a configuration. We find that perturbations 
resulting from the quenches propagate at smaller group velocities as the gap controlling the spatial extent of the 
edge modes decreases. Our results provide different routes to check dynamically the non-trivial nature of Chern 
bands. 


I. INTRODUCTION 


Chern insulators are topological states of matter that host an intrinsic quantized quantum Hall effect without the need of 
external magnetic fields. Embedded today in the large body of knowledge gathered around topological phases IHIIl, Chern 
insulators were conceived soon after the discovery of the quantum Hall effect a and can hold the key to realizing interesting 
unconventional phenomena such as fractional many-body states in the absence of external magnetic fields ilia or more exotic 
electrodynamical phenomena such as a repulsive Casimir effect a. 

The equilibrium signatures of Chern insulators are by now well established. Their anomalous transport is determined by a 
topological invariant associated to each band known as the Chern number related to the Hall conductivity through Q axy = 
Cae^/h, where Ca is the Chern number for each filled band a. The non-trivial band topology, i.e. non-zero Chern numbers, 
implies the existence of propagating chiral modes localized at the edge of a finite sample. The recent experimental realization of 
solid-state materials iHni and optical-lattice systems GMll having Chern bands (or closely related properties) has boosted 
the interest in such systems. 

Compared to equilibrium behaviors, out-of-equilibrium features of Chern bands have attracted interest only recently, partic¬ 
ularly because realizations with ultracold atoms are well-suited for studying real-time dynamics. It is therefore experimentally 
relevant to characterize topological states of matter not only with transport signatures but also with dynamical measurements ESl- 
1^ . For instance, the Chern number can be determined dynamically by following the trajectory of a wave packet in a tilted optical 
lattice, MM a method which was successfully used experimentally recently. iniiiii It was also shown theoretically that edge 
state dynamics can be visualized by tailoring external potentials that, when released, provide direct imaging of the chirality of 
the Chern bands through the chiral propagation of the atoms in a hard-wall trap ll20l . Global quenches in Chern band models 
and related systems have been addressed in a number of recent studies ll2Tl - E9l . 

In this work, we complement this body of knowledge by studying the effect of local quenches on Chern bands filled with 
fermions. We focus particularly on the real-time chiral response of the edges. This study thus supplements existing proposals 
for experimental protocols to probe the chiral nature of such states ll23l - E5l 130143^ . 

We consider two kinds of quantum quenches shown schematically in Fig. a single quench and a pulse consisting of two 
sudden changes of the Hamiltonian. The background Hamiltonian Hq is taken to be a tight-binding model with topologically 
non-trivial bands. The perturbing Hamiltonian Hi is chosen to be a localized density perturbation of the system. We will 
concentrate on perturbations at one site at the edge of the fermionic system, in order to dynamically probe the chirality of the 



FIG. I. (color online) The figures schematically show the two quench protocols considered in this work: (a) is a quench from an equilibrium 
state att = to while (b) is a pulse of width t = t2 — t\. 














2 



FIG. 2. (color online) The single-particle energy spectrum for the model 0 on a 17 X 17 lattice (each site has two orbitals) with open boundary 
conditions. The parameters are chosen to be in (a) the Chern insulator phase and (b) the trivial insulator phase. The insets show the real space 
probability distribution of the single-particle eigenstates closest to half-filling (marked by the arrows on the spectrum) with the color intensity 
showing the probability at any particular site. 


edge states. In the first case, we start with the ground state of and then remove the local perturbation Hi and follow time 

evolution under the Hamiltonian Hq. In the second case we start with the ground state of Hq and turn on the local perturbation 
Hi for a finite time interval r. 

We first consider (section[In| half-filling of an open-boundary system. In this case, the edge of the fermionic system is simply 
the edge of the lattice. Next, motivated by experiments with cold atomic systems, we explore the effect of power-law traps 
(section [Tvl, i.e., traps of the form ~ . In cold-atom experiments, the atoms are usually loaded in a harmonic trap, i.e, 7 = 2 . 

Since systems in harmonic traps often do not have a well-defined sharp edge, there is significant interest in power-law traps with 
large exponents 7, which is expected to be more similar to a system with sharp boundaries lEollIIlETlEsl. Accordingly, we 
consider fermions in a Chem band in the presence of a power-law trapping potential J(r/ro)''', where J is the hopping energy 
that sets the energy width of the bands. The parameter tq then functions as a trapping length scale, determining the spatial extent 
of the trapped fermionic cloud lITSl 1^ |2T1 . We first characterize the equilibrium properties in such a trap and then study the 
effect of local quenches on it. 


Our main findings are as follows. When the bands are topologically non-trivial, a local quench of a site potential at the edge 
generates a spatially localized pulse that propagates chirally around the edge of the sample. This is a straightforward and very 
direct manifestation of the topological nature of Hq. The front of the pulse has speed determined by the hopping scale J, but 
the peak of the pulse has smaller speed for smaller bulk gaps at the relevant points in the Brillouin Zone. We attribute this 
phenomena to the increasingly poor localization of the edge states as the gap decreases; the perturbation is spread out over 
increasingly more lattice sites rendering the sharp boundary picture less accurate. In the second part, considering power-law 
traps with large exponents 7, we show how the trap parameter tq determines an optimal particle number for mimicking the 
half-filling situation, and hence for observing the chiral nature of the edge states. 


The paper is structured as follows. In sectionjl^we present a generic model defined on the square lattice that has both trivial 
insulator and Chem insulator phases; we use this model throughout the paper. Next, in section [1111 we consider a half-filled 


rectangular lattice with open boundaries and study local perturbations at the edge, of the two forms shown in Fig. [T] (quench 
and pulse). In section|IV]we consider power-law potential traps. We exemplify the decisive role of rg in selecting the optimal 
fermion number in order to mimic a half-filled system with a populated edge. We show how the non-equilibrium features (chiral 
propagation) survive in this case. In section |V] we summarize and provide context, and discuss the experimental set-ups that 
could access the chiral nature of the Chern bands as proposed here. Since the Hamiltonians are non-interacting, the methods for 
finding ground states and time-evolving many-body states are standard, and are briefly described in Appendix [A| 


II. MODEL HAMILTONIAN AND QUENCH PROTOCOLS 


We are interested in the dynamics of a generic two band Chern insulator model in two dimensions. As a representative of 
a Chem-band system, we will consider the following single particle Hamiltonian defined on the square lattice ll39l l40l with a 













spin-1/2 degree of freedom (either orbital or spin) 

Hci = Hj + Hm, 
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Hm = Mc\azCi, (lb) 


with hopping J between the sub-lattices and a staggered chemical potential M. The fermionic operators c^, for each site i 
represents a two-spinor of the operators for each orbital/spin on that site, e.g., c\ = . Summations over the spin index 

are thus implicit. 

In Fourier space, a two band model Hamiltonian can be written generally as 


iTci(k) = ek -f cr • dk. 


( 2 ) 


In our case, dk and Ck 

= - J 

dy k - J Sin(/Ly), 

rfz,k = M — J ^ cos{ki), and Ck = 0. 

i—x,y 


(3a) 

(3b) 

(3c) 


At half-filling such a single particle Hamiltonian represents a Chern insulator with chiral edge states whenever —2 J < M < 2J. 
The topological and trivial eigenstates for a finite sample are shown in Fig.|^(a) and (b) respectively, together with the real 
space probability distribution of the eigenstates closest to zero energy {E = 0). In the trivial case, there is no sign of edge states 
whereas for the Chern insulator, the edge states are clearly visible. At half-filling the Chem insulator state has axy = Ce^ jh 
where C = ±1 is the Chern number of the filled lower band. We note that Ck = 0 does not affect the topological nature of the 
state and thus it is irrelevant for our discussion. 

The Hamiltonian Q (with small modifications) has been used widely as a basic standard model for two-band systems with 
each model having a Chern number Il39l - l4n . The Haldane honeycomb model lO can be cast into this form if the two sites in 
the honeycomb unit cell are mapped into the two orbitals on the same site of Hamiltonian Q. The results we present are quite 
generic, and should hold qualitatively for other Chern insulator models such as the Haldane model 0, and also for other systems 
having bands with nonzero Chern number, such as the Hofstadter model 1421 . 

In section III we will set Hq = Hci while in section IV a trapping potential iTtrap (to be defined later on) will be added to the 
HamiltonianrHo = Hci + Tftrap- 

We will be concerned with physical processes where the system described by Hamiltonian Hq is acted upon by two types of 
perturbation. The first is of the form 


H,{t) = H,9{to-t), (4) 

which we will refer to as a quench from equilibrium [see Fig.[2a)] and the second is defined as 

Heit) = Hi[eit-ti)-0it-t2)], (5) 

that is non-zero for a time t = t 2 — ti, and we refer to as a pulse [see Fig.[2b)]. We focus on local density perturbations that 
we label Hi. We define a local density perturbation as an increase or deficit of the charge density around a particular site I with 
magnitude p/ 


H[=p,ic]ci, (6) 

Here no implicit summation is assumed over the site index, but there is an automatic summation over the spin index, since c; 
and c| are two-spinors. 

In the concluding section we will discuss briefly possible physical implementations of these types of perturbations in cold- 
atomic experiments and in condensed matter settings. 

Next, we perform time evolution by evolving the one-particle density matrix, or the matrix of correlators, Qij = (c|cj) (see 
Appendix A for details); here i, j are site indices and the spin indices are implicit. We present our results by plotting the total 
densities at each site i, which correspond to the diagonal terms of the correlation matrix, g^, relative to gufi, the density at site 
i calculated from Hq, with spin summation implied. 
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FIG. 3. (color online) The time evolution of the density in the Chera insulator, following a quench from equilibrium [Fig.[TJa)] with to = 0. 
Snapshots of the local local density difference gait) — gufi at each site i are shown at four different times. A red (blue) color on any site 
denotes an excess (deficit) in fermion density relative to ga^ with the intensity of the color showing the magnitude of the difference. The 
arrows show the direction of the chiral propagation of the excitations. The local quench is performed at the site labeled by I and the density 
measurements in Fig. |^are done at the site labeled n. For this simulation/r; = —J/4and M/J = 1. Please see supplementary material for a 
video of the time evolution. 


III. HALF-FILLED OPEN-BOUNDARY SYSTEM 


In this section we consider rectangular open-boundary systems (without a trap) with L 
explore the chirality of the topological edge states through perturbations of the form 


Lx X Ly sites at half filling. We 


A. Local density quench from equilibrium 

We start with a local density quench, of the type of Fig. [TJa), at an edge site I, 


H{t)=Ho + H[e{to-t), (7) 

with Hq — Hqi and H\ given by (0. Through the action of H[ the initial ground state has different density at site I G edge 
compared to the ground state of Hq. The perturbation p/ has equal effect on both spin components. The particle number is fixed 
to N = Lx X Ly, i.e. half-filling. 

In Fig.|^we show the time evolution by presenting four snapshots of the local local density difference gait) — qh^q between the 
density at each site i (gii{t)) and the one corresponding to the half-filled ground state of Hq (ga^o). It is clear from the snapshots 
that the perturbations travel along the edge of simulated lattice with a well-defined chirality, in this case counterclockwise. The 
effect of the perturbation on the edge site labeled n, marked in Fig. is shown in Fig. The density difference at site n shows 
a periodic pattern with a period Tedge — P/vp, where P is the perimeter of the lattice and uf is the Fermi velocity of the edge 
state. The perimeter for this particular simulation is P = 16a x 4 = 64a. 

In order to shed light on these numerical results, we estimate the value of vp and the localization length of the edge from 
the low energy theory of the model ([T]|. To do so, we first expand Hq around each high-symmetry point of the Brillouin Zone, 
k(* i) = TT{i,j) with i,j = 0,1. The emergent effective low energy model is a massive Dirac equation around each one of them 
that takes the form 


= -J [(-l)*fc, + {-lyky] + 


( 8 ) 
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FIG. 4. (color online) The density difference Qnn{t) — Qnn,o at the site n marked in Fig.j^is plotted against time for a trivial insulator with 
M/J = 3 (featureless green solid line), a Chern insulator with M/J = 1 (blue solid line) and a Chem insulator with M/J — 1/2 (orange 
dashed line). The gap sizes at k = 0 are respectively 2\M — 2J\ = 2J, 2J and J, the local density perturbation parameter is /r; = — J/4. 
Note the slight delay for the dashed curve, corresponding to a smaller gap at k = 0 and hence more extended edge state (see Eq. l[^). The 
time interval Tedge — P/vf, where P = 64a is the perimeter of the lattice and vf = Ja/h is the Fermi velocity of the edge state, is also 
shown by the black dashed lines. 


with four mass terms defined by 

= M - [i-iy + i-iy] J. (9) 

If \M/ J\ < 2, i.e. in the topological phase, only three out of the four mass terms have the same sign. On the other hand when 
\M/ J\ > 2, the trivial state, all of them have the same sign. For M/ J > 0 (< 0) the boundary between the Chern insulator and 
a trivial insulator (e.g. vacuum) is modelled by choosing TO(o_o) = tn(t/) = rn{y)) such that m{y) changes sign at the 

boundary, which we take it to be at y = 0. The corresponding Dirac equation has only as a good quantum number and has a 
gapless solution that decays exponentially as ll43l 


'I'edge(2/) 


g- /o” m.{y')dv' 


1 

1 


( 10 ) 


Fixing M/J > 0, a sharp edge can be modeled by m(y) = [9{y) — 6{—y)], that determines the localization of the edge 

state to be inversely proportional to \M — 2J\ (the case where M/J < 0 is obtained by simply replacing m(o,o) ”^( 7 r, 7 r))- 
Such an edge state, has a dispersion E = Evpky (E = Evpkx) for edges along the y (x) direction. The sign is determined by 
the sign of the Chern number C of the lower band and wj? is set by the bulk dispersion Fermi velocity. Therefore, vp can be read 
directly from ( [Saj i. For this model it is isotropic and takes the value vp = J (in units of K). 

From this analysis it follows that, if is reduced, the edge states will have a hnite extent, having in general support on 

several rows close to the edge. This in turn will affect the period between the density pulses reaching a particular site, i.e., the 
period between the peaks in Fig. Since the edge perturbation now has more sites to explore as it propagates, one expects that 
the front of the propagating wave travels at the same speed as in a narrow edge, but the peak of the density wave will travel 
more slowly due to the larger width of the propagation channel. This effect is shown in Fig. where the local density at a site 
n € edge given by gnn(t) is shown for two different instances within the Chern insulator phase, corresponding to M/J = 1 
(blue-solid line) and M/J = 1/2 (orange-dashed line). Fig.|^shows how the highest crest of the oscillations shift to later times 
as the \M — 2 J| decreases from J to J/2 as argued above. The speed of the front of the wave is apparently unaltered as expected 
because vpis independent of TO(o,o)- 

Such an effect is observable in a simple toy model of the conducting chiral edge by analyzing how a perturbation propagates 
along a conducting strip modeled as a trivial tight binding square lattice of linear dimensions x ly with ^ ly, as a function 
of the strip thickness ly. Although in this case the propagation is not chiral and thus an exact comparison is not possible, the 
wave packet indeed explores more sites as ly is increased, which effectively reduces the peak velocity in the x direction, similar 
to what we observe for the Chem insulator edge. 

The apparent decay of the main peak in Fig.j^on the other hand is likely to be due to the spread of the wave packet along the 
edge rather than decay into the bulk. The latter is strongly suppressed by the initial perturbation being localized at a single edge 
site, having therefore no overlap with bulk states. 
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FIG. 5. (color online) The local density difference Quit) — Qufi is plotted for different times to show the dynamics of the edge excitation in 
the case of the pulse quench [Fig[^b)] with fi = 0 and t 2 = 30h/ J. The red and blue arrows show the chiral propagation of the particle- and 
hole-like excitations respectively. The quench protocol parameters are chosen such that r > h/ J. Other parameters are the same as in Fig.j^ 
Please see supplementary material for a video of the time evolution. 


Fig. I^also shows data for the trivial insulator case, M/ J = 3 (featureless solid green line). There is no chiral propagation in 
this case due to the absence of edge states, so there are no features observed in the density at site n. 

Note in addition that in Figs.j^and Hthe perturbation satisfies < \M — 2J\. We have checked as well that our results hold 
for different values of When fxi >\M — 2J\ the oscillations are more pronounced due to the enhanced finite overlap of the 
initial perturbation with bulk states. Otherwise, the conclusions stay as above. 

Finally, we have also studied the decay of a local perturbation at a bulk site, as a function of the gap size. The resulting 
disturbance spreads out in all directions and there is no notion of chirality in the propagation of the excitation. We also find 
no correlation with the gap. We attribute this to the fact that a spatially localized bulk perturbation is strongly delocalized in 
momentum space and thus it is insensitive to the size of the gap. 


B. Local density pulse 

We now study a density perturbation pulse [see Fig. [TJb)] of width r, 

H{t)=Ho + H[[e{t-t,)-9it-t2)], (11) 

with Hq = Hci and t = t 2 — ti. We choose ti = 0 and t 2 = t without loss of generality. 

We focus on a local edge perturbation (at a site I G edge) and look at the time evolution as a function of r. Since the relevant 
energy scale is J we expect two different regimes depending on whether r < h/ J or r > h/ J. Indeed, in the former case 
(r < h/ J) the perturbation is too localized in time to react separately to both the turning on of the perturbation at ti and 
releasing at ( 2 - Accordingly, the system generates a single signal traveling along the edge of the system, which is allowed to 
conduct chirally due to the non-trivial topology of the Hamiltonian. Such a disturbance evolves in time qualitatively as was 
shown in Fig.j^and thus is not shown here. 

On the other hand, when t > H/J the system can react to both the pressing at ti and the release at ^ 2 - In this case, two 
pulses are generated, one at time ti when the system is ‘pressed’ (perturbation is turned on) and one at (2 when the local density 
perturbation is released [see Fig.[^. 

The Fermi velocities of both pulses are set by uf = J since the arguments presented in the previous section still apply. 
Similarly, by changing M [see Fig.|^ a slight delay in the wave packet center propagation, analogous to that observed in Fig.|^ 
is evident. 
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FIG. 6. (color online) The density difference — g„„,o at the site labeled n [see Fig.[^ is plotted against time for the case of the pulse 

quench with = 0 and t 2 = 30h/ J for parameters M/J = 3 (featureless green solid line), M/J = 1 (blue solid line) and M/J — 1/2 
(orange dashed line). As in Fig. 0 there is a slight period delay as the gap decreases due to the finite extent of the edge state. The intervals 
T = t2 — ti as well as Tedge — PJvf, where P is the perimeter of the lattice and vp is the Fermi velocity of the edge state, are indicated by 
vertical dashed lines. 


IV. CHIRAL ‘EDGE’ DYNAMICS IN POWER-LAW TRAPS 
In this section, we consider fermions in a Chern-band lattice in the presence of a power-law trap 


Ho = Hci + i7trap, ( 12 a) 

Htr.p = jY.(^) (12b) 

i V ^0 / 

where is the distance of site i from the center of the trap, which in our simulations will coincide with the center of the lattice. 
The purpose of considering a trapped Chern lattice is to make contact to their possible cold-atom realizations. 

Cold atom experiments are generally performed in the presence of harmonic traps, i.e., 7 = 2 traps. In such a trap, it may be 
difficult to distinguish between edge and bulk regions. As a result there is considerable interest in power-law traps with larger 
exponents, i.e., large 7 , because for large 7 a power-law trap can resemble a box trap with sharp boundaries, and hence can 
be expected to have better-defined edge states spatially separated from the bulk ll20l 1271 [37l l3^ . For example. Refs. l20l and 
iTTInroDose non-equilibrium protocols involving hard boundaries for experimentally determining Chern numbers, using the idea 
that approximate hard-wall potentials can be designed using power-law potentials with large 7 . 

Motivated by this interest in large 7 , in this section we present equilibrium spatial structures of fermions loaded in a 7 = 50 
trap, and based on this knowledge we explore chiral responses of the ‘edges’ of such fermionic clouds to local quenches. 

The effect of a large -7 trap is to separate energetically single-particle eigenstates which are spatially ‘inside’ the trap (within 
distance tq from the trap center) from those ‘outside’ the trap. A low-energy sector thus emerges which closely resembles 
a uniform system of radius tq. By ‘half-filling’ this region, i.e., by having N ~ TTr^/a^ fermions in the entire system, one 
can then mimic a half-filled region with a reasonably well-defined edge. Such an ‘edge’ also shows chiral dynamics as in the 
open-boundary case without a trap treated in the previous section. 

In Fig.[^ we use a square lattice with sides larger than 2ro, and 7 = 50. The top panels show the single-particle eigenspectrum. 
The higher energy sectors contain eigenstates whose weights are spatially concentrated in regions r > rp. The zoom onto the 
low-energy sector [see top right panel in Fig. [7] resembles the spectrum of a Chern lattice without a trap and with open boundary 
conditions. Indeed, eigenstates with —I<i7/J<l have an edge-like distribution localized around r ^ tq, exemplified by the 
lower right panel in Fig.|^ On the other hand, eigenstates with 1 < \E/J\ < 3 are extended over r < ro, resembling the bulk 
behavior of a Chern insulator [see bottom left panel in Fig.|^. 

Note that the midgap ’edge’ states with —1 < E/J < lare not as sharply defined as in the free Chern insulator case, as 
can be seen by comparing the lower panels in Fig. [^with Fig. |^(a). This is because, even at such large 7 , the lack of a hard 
wall induces mixing between eigenstates with edge and bulk character. Nevertheless, the edge modes are well-defined enough 
to display chirality in real-time dynamics, as we next show. 
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FIG. 7. (color online) Single-particle energy spectrum and eigenstates in the presence of a trap, ro = 6 a and 7 = 50, for M/ J — 1 
and Lx ~ Ly = 17. (a) Energy spectrum. Inset shows absolute values on a logarithmic scale, (b) Magnification of low energy sector. 
The eigenvalues in this sector resemble that of a Chem insulator with open boundary conditions [compare with Fig.|^a)]. (c,d) Probability 
distributions for the a = 60 and a = 120 eigenstates (dashed lines in (b)), corresponding to a ‘bulk state’ and an ‘edge state’. 



FIG. 8 . (color online) The density along the x direction for a fixed cut y = 8a in real space for different values of ro for the ground state with 
N — 120 ~ 7 r 6 ^ fermions, 7 = 50, M/J — 1 and Lx = Ly = 17. For this N, the edge states are only visible for ro = 6 a as bumps in the 
density at the edge of the trapped cloud. 


To access the chiral character of these effective edge states, it is necessary to mimic the situation of half-filling. This is achieved 
by having the fermion number to be close to . In Fig. |^we fix the fermion number to be iV = 120 ^ ttG^ ~ 113, and 

compare the density profiles along the horizontal (x—) direction for traps with ro = 5a, ro = 6a, and ro = 7a. The figure shows 
that this particular filling provides enough fermions to occupy the effective edge states when ro = 6a. The edge occupancy 
is only visible in the ro = 6a case through the bumps in the density profile peaked at x = 3a and x = 13a; the other curves 
lack this feature. This exemplifies the fact that, depending on the trap shape (set by ro), there is an optimal particle number for 
accessing the chiral edge states of the system. 
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FIG. 9. (color online) Four snapshots of the density perturbation Qii{t) — gufi after a quench of the form of Fig.[^a), for the Hamiltonian Hci 
with \M\ < 2 J in the presence of a 7 = 50 trap. Trap radius is tq = 6a (upper panels) or ro = 7a (lower panels). The simulation lattice is 
17a X 17a; the quench strength is = — J/4. The particle number N = 120 « is chosen such that the last filled eigenvalue for ro = 6 
is that of Fig. |^d). Here ga^ corresponds to the site densities in the ground state of Hq = Hci + Htrap filled with = 120 fermions. Please 
see supplementary material for a video of the time evolution. 


In order to highlight the importance of the filling fraction with respect to the trap parameters, in Fig.|^ we show the real-time 
evolution to a perturbation at a fixed site I using the non-equilibrium protocol of Fig.[TJa) for tq = 6 a and tq = 7a. The particle 
number is fixed for both cases to be (V = 120 ^ ttG^ « 113, the same value as in Fig. Therefore the edge states are only 
populated for tq = 6 a but not for rg = 7a. For rg = 6 a (upper panels of Fig. |^, the perturbation follows the profile of the 
trap in a chiral fashion, a situation similar to a half-filled open-boundary system. For rg = 7a, the density disturbance does not 
propagate along the new edge, but dissipates into the bulk. 

Note that there is some flexibility in varying N for a fixed rg (and vice versa) while still obtaining a visible edge state. As 
long as the state at the Fermi energy is one of the edge states, the physics of edge states is accessible. The number of edge states 
is approximately 27rrg, the circumference of the interior region. The lattice geometry will affect the exact number, of course, 
but for rg ^ 1, this is a reasonable estimate. For smaller rg, the spectrum in Figure 7 exemplifies the situation: in this case 
(rg = 6 a), the edge states run from eigenstate a ~ 95 to eigenstate a ^ 131. The physics of edge states is visible as long as N 
is in this range. 

For the rg = 6a case, as the gap is decreased by changing the ratio M/J appropriately, the effective edge states become less 
localized. Fixing the rest of parameters we observe a retardation effect of the wave packet, consistent with that discussed in Sec. 


Ill for the evolution in the absence of Ht 


trap* 


For completeness we have also investigated the protocol in Fig. 1 (b) under the effect of the trap and found similar results 
to the scenario without a trap as long as the system is close to the optimal Ailing discussed above. This protocol generates two 
pulses confined to the boundary of the trap with similar properties as those discussed in previous sections. Finally we expect that 
the effect of softer traps, i.e. smaller values of 7, addressed for instance in Refs. |20] andwill result in an overall broadening 
of the edge states and greater mixing between bulk and edge. To what extent the chirality in real-time dynamics is visible for 
smaller 7 remains an open question. 


V. DISCUSSION AND CONCLUSIONS 


In this work we have investigated the effect of local edge quenches on the dynamics of Chern bands. We have addressed two 
types of local density quench protocols to probe directly the chirality of the edge states, both with and without confining traps. 
When the bands are topologically non-trivial, these quenches generate localized edge excitations that propagate chirally along 
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the sample with a group velocity that decreases as the gap controlling the spatial delocalization of the edge modes decreases. In 
addition, we have shown that, in the presence of a sharp trap with a characteristic length scale rp, the chiral dynamics is only 
observable if the filling is such that the effective trap confinement region is half-filled. These results emphasize the importance 
of the effective filling within the trap to observe edge dynamics, even when the trap potential is quite sharp. 

We have concentrated on large-exponent power-law traps (large 7). While this is not the common situation in current cold- 
atom experiments, interest in producing and utilizing such traps is high, particularly in the context of topological matter ll20l 
IZU |32l m. We have shown data for 7 = 50, but it is expected that traps with exponents down to 7 « 4 will have similar 
properties ll44ll . By analyzing the spectrum, we have shown that the ‘inside’ and ‘outside’ of the trap are energetically separated. 
The ‘inside’ region is found to cover Trrg lattice sites. This observation has allowed us to specify the particle number which 
mimics the physics of half filling in the ‘inside’ region. We have shown that the chiral dynamics is not visible when the number 
of fermions is very different from this optimal, because the effective filling then differs from half-filling. 

It is plausible that the protocols considered in this work can be realized in cold atomic set-ups. Addressing and imaging 
these systems with single-site resolution have become available in cold-atom laboratories during the past few years msizi. 
Thus, performing local quenches on the recently realized Chern-band lattices ifTTlfTSlISTl and following the ensuing site density 
dynamics in real time should be technologically feasible. 

Our findings might also be relevant in solid state set-ups. For instance, the Chern insulator state has been recently realized by 
magnetic doping a thin-film structure of a three-dimensional topological insulator m. In this context, the protocol in Fig.[2a), 
theoretically can be thought of as an STM tip perturbing the system locally to then release the density. The typical time scales 
governing such dynamics in this case are several orders of magnitude faster than in cold atomic experiments, rendering such a 
proposal practically unviable experimentally. However, fast, out-of-equilibrium photo excitation of electrons has been measured 
with recent pump-probe techniques HSi . In this experiment, the local dynamics occurring within typical electron time-scales 
of femto-seconds were measured at the surface of a three-dimensional topological insulator. In light of these results we can 
reinterpret the protocol in Fig. [TJb) in a first approximation as a laser pulse of duration r that probes the electrons locally in 
space and time. It is therefore not unrealistic to think that implementing such a protocol might be possible by probing the Chern 
insulator state of Ref. [H especially considering that they are grown from essentially the same family of materials that were 
pump-probed l48l . 
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Appendix A: Density matrix time evolution 

Time evolution in this work has been performed by evolving the single-particle density matrix, i.e., the matrix of two-point 
correlators. The method is standard for systems described by quadratic fermions and is widely used in the non-equilibrium 
literature, so we outline the procedure here only briefly. The initial equilibrium density operator is p = IV'i) (V'il’ where the 
sum is over occupied single particle eigenstates \tjji) of the initial Hamiltonian. 

The time evolution of such an operator is given by 


g(t) = (Al) 

where H is the single-particle Hamiltonian under which the evolution occurs. As single-particle operators, the density and 
Hamiltonian operators can each be represented as 2 x x Ly matrices in our case. For the lattice sizes we have used, the 
matrix dimensions are < 10^; the matrix operations of Eq,. ( |A1| | are thus numerically inexpensive. The g{t) matrices at each 
time can thus be performed either by rotating to the basis where H is diagonal, or even by explicit computation of the matrix 
exponentials. For the pulse case, the Hamiltonian changes at the end of the pulse; we perform evolution with Hamiltonian 
Hq -f Hi for time r and then switch to evolving with the Hamiltonian Hq. 

Physically, the diagonal elements of gui^t) represent the evolution of the density at site i with the constraint that Tr[p(f)] = N 
where N is the number of particles. The off diagonal elements gij(t) with i j are equal-time correlation functions at time t 
between sites {i,j}. 
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